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Abstract 


The  distortion  of  the  growth  rates  of  unstable  baroclinic  waves  is 
investigated  using  implicit  finite  differences  to  replace  the  linearized 
governing  equations  of  a  two-parameter  model.  A  variety  of  space  and  time 
mesh  sizes  Ax  and  At  are  used,  together  with  three  values  of  the  static 
stability  parameter  d&/dp.  On  the  whole,  for  commonly  used  values  of 
Ax  and  At,  and  for  typical  values  of  d0/3p,  the  introduction  of  finite 
differences  appears  to  cause  only  a  slight  distortion  of  the  growth  rates 
over  the  major  portion  of  the  unstable  spectrum,  the  long  waves  being  less 
affected  than  the  short  waves.  Effects  due  to  At  are  particularly  minor. 
The  finite-difference  growth  rates  depend  very  slightly  on  the  zonal  wind 
U,  in  contrast  to  the  continuous  case  where  no  such  dependence  is  found. 


1,  Introduction 


In  a  recent  paper  Gates  (1961)  discussed  the  effect  of  finite- 
differences  on  baroclinic  instability  criteria  for  the  case  of  a  linearized 
two-parameter  model,  using  an  implicit  difference  scheme.  His  results  show 
that,  for  a  given  vertical  windshear,  the  spectrum  of  unstable  waves  is 
shifted  slightly  toward  shorter  wavelengths,  leaving  the  minimum  critical 
shear  unaltered.  The  magnitude  of  the  shift  was  found  to  be  directly  pro¬ 
portional  to  the  space  mesh,  Ax,  and  inversely  proportional  to  the  static- 
stability.  In  this  connection  it  might  be  interesting  to  know  what  effect, 
if  any,  the  introduction  of  certain  finite  differences  has  on  the  growth 
rates  of  unstable  baroclinic  waves.  It  is  the  purpose  of  the  present  dis¬ 
cussion  to  throw  some  light  on  this  problem.  The  same  atmospheric  model, 
and  the  same  differencing  scheme  as  in  the  paper  mentioned  above,  will 
be  used. 


2.  Baroclinic  difference  equations 


The  particular  baroclinic  model  employed  is  the  thermotropic  model, 
whose  details  are  described  elsewhere  (Thompson  and  Gates,  1956),  and  whose 
linearized  governing  differential  equations,  with  the  usual  assumption  of 
y-independent  flow,  may  be  written  as 
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where  f  and  4  are  the  streamf unctions  of  the  vertically-integrated  flow 
and  the  thermal  wind,  respectively.  U  is  the  zonal  wind,  and  U*  is  the 
shear  or  thermal  wind,  both  assumed  constant;  B  is  the  Rossby  parameter, 
and  depends  on  the  static  stability,  and  is  given  by 
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where  f  is  the  Coriolis  parameter,  R  the  specific  gas  constant  for  dry  air, 
p0  *  1000  mb,  T  is  the  absolute  temperature,  and  0  is  the  potential  tem¬ 
perature.  For  the  purpose  of  avoiding  extraneous  roots  in  the  frequency 
equation  below,  and  of  using  a  numerically  stable  scheme,  implicit  differences 
are  now  introduced,  and  the  appropriate  finite-difference  analogues  of  (1) 
and  (2)  are 
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where  t  3  t(mAx, TAt) ,  i  _  =  $(mAx,TAt),  and  where  Ax  and  At  denote 
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space  and  time  increments,  respectively. 

In  order  to  investigate  the  behavior  of  growth  rates,  it  is  con¬ 
venient  to  assume  that  (3)  and  (4)  have  solutions  of  the  form 


where  i  is  the  imaginary  unit,  a  *  kAx,  n  =  kC^,  k  *  2*/L,  and  where  L  is 
the  wavelength,  and  is  the  phase  speed  of  the  finite-difference  solution. 
The  amplitudes^  and  £  are  assumed  constant.  Substituting  (5)  into  (3) 
and  (4),  one  finds  readily  enough  that  the  system  of  equations  for  the 
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un known  amplitudes  is 
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The  condition  that  (6)  should  have  non-trivial  solutions  requires  that  the 
determinant  of  the  coefficients  of  £  and  £  must  vanish.  This  leads  to  the 
frequency  equation  of  the  numerical  solution f 
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In  comparison,  one  finds  by  the  usual  methods  that  the  frequency 
equation  for  the  differential  system  (1)  and  (2)  is 
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where  C  is  the  phase  speed,  and  k  =  2n/L  as  above.  Equation  (9)  expresses 
the  well-known  fact  that  the  numerical  phase  speed  Cu  depends  on  the  choice 

I* 

of  Ax  and  At. 


3.  Growth  rates  in  the  difference  case 


Let  CN  be  complex  with  real  part  C„  ,  and  imaginary  part  CN  ,  i.e., 

r  i 

CM  *  CM  +  iCM  .  Since  (5)  is  a  solution  of  the  difference  system  (3)  and 
r  i 

(4),  subject  to  condition  (9) ,  unstable  baroclinic  waves  can  exist  in  this 
system  if  Cj,  has  a  non-trivial  imaginary  part.  If  (5)  is  now  written  in 
the  form 


(12) 


it  is  clear  that  the  waves  will  amplify  provided  CN  >0,  Now,  the  inverse 

l 

tangent  in  (9)  has  real  solutions  for  all  real  values  of  XZf  and  can  be 

complex  only  if  Z  is  complex,  because  X  is  always  real.  For  complex  Z, 

let  Z  a  Z  +  iZ.  where,  from  (10) 
r  x 
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Thus,  the  critical  curve,  that  is,  the  curve  which  separates  the  region 
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of  stability  from  that  of  instability,  is  defined  by  the  requirement 
Z^  =  0.  The  analogous  curve  in  the  differential  system  is  given  by  the 
condition 

-  (■£  -ph?)  - 


which  is  the  same  as  setting  Z^  =  0,  and  dividing  xa  by  The 

cut-off  wavelength  L  ,  such  that  all  waves  for  which  L  £  L  are  stable, 

c  c 

the  minimum  shear  U*  such  that  all  waves  for  which  U*  £  U*  are  stable, 
mf  o  f 

and  the  wavelength  L  corresponding  to  U*  are  all  obtainable  from  (14), 

m  m 

Expressions  for  L  ,  U* ,  and  L  have  already  been  given  elsewhere  (Gates,  1961) 
c  m  m 

and  it  need  only  be  mentioned  here  that  the  method  used  in  the  present 
paper  gives  the  same  results  for  these  parameters* 

To  discuss  the  effect  of  finite  differences  on  the  growth  rate,  it 
is  necessary  to  find  the  imaginary  part  of  the  phase  speed  (9).  Thus, 
letting  Z  be  complex,  one  finds  after  some  manipulation  that  (See  Appendix), 
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The  form  of  (16)  indicates  that  the  positive  square  root  should  be  taken 

in  (14)  for  amplifying  waves  (C  >  0),  Similarly,  CN  <  0,  and  the  waves 

Nt  wl 

decay  if  the  negative  root  is  taken  in  (14),  For  finite,  non-zero  Ax  and 


-7- 


At  the  waves  are  neutral  (C^  =  0)  when  Z^  =  0,  and  when  X  =  0.  The  case 
Zj  =  0  is  the  one  mentioned  above.  The  case  X  =  0  occurs  whenever  Ax  is  any 
integer  multiple  of  L/2.  This  means,  at  least  in  the  linear  case,  that 
the  harmonic  components  whose  wavelengths  are  below  the  resolution  limit 
of  the  grid  remain  reutral  when  an  implicit  difference  scheme  is  used. 

Various  expressions  for  the  growth  rate  are  possible;  but  in  the  finite- 
difference  case,  a  convenient  measure  of  this  quantity  is  the  doubling  time. 
This  is  the  time  required  by  the  wave  amplitude  to  double  in  magnitude, 
and  can  be  written  as 


-A2 


(17) 


where  CN  is  given  by  (16).  The  distortion  of  the  doubling  time  by  implicit 
finite  differences  is  shown  in  figs.  1-9  for  selected  values  of  Ax,  At, 
and  the  static  stability  0^dp.  It  is  evident  from  the  figures  that  the 
distortion  is  slight  for  commonly  used  values  of  Ax  and  At,  and  for  a 
typical  value  of  dE/dp  (-6  deg/100  mb).  These  figures  also  show  that 
the  distortion  is  more  serious  for  large  values  of  Ax,  short  waves,  and 
small  static  stabilities,  in  agreement  with  the  behavior  of  the  shift  of 
baroclinic  instability  criteria  found  by  Gates  (1961).  In  general,  the 
growth  rate  of  the  unstable  shorter  waves  is  increased,  and  that  of  the 
longer  waves  decreased,  except  in  the  case  of  neutral  static  stability 
(dfl/dp  *  0) ,  where  the  growth  rate  is  always  decreased.  Another  effect 
of  the  finite  differences  is  that,  for  a  given  shear  U*,  the  wavelength 
of  the  most  rapidly  growing  wave  is  slightly  lowered.  This  shift  is  more 
pronounced  for  the  smaller  shears,  larger  space  mesh  sizes,  and  small 
static  stabilities.  For  a  given  static  stability,  the  maximum  value  of 
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Figure  1,  Doubling  time  in  days.  The  solid  lines  are  for  the  finite- 
difference  case,  the  dashed  lines  for  the  continuous  case.  The  latitude 
is  45  deg  N.  Here  the  static  stability  parameter  d6/3p  =  0,  Ax  =  3OO  km, 
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e  2.  Same  as  fig.  lf  but  with  Ax  -  6 00  km. 
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Figure  3.  The  effect  of  the  space  mesh  Ax  on  a  doubling  time  of  exactly 
1  day.  The  latitude  is  45  deg  N,  dfc/dp  =  0,  and  At  =  1  hr. 
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this  shift  is  equal  to  the  displacement  of  the  minimum  point  on  the  critical 
curve,  corresponding  to  a  given  Ax,  and  decreases  with  increasing  shear. 
Typical  maximum  shifts  are  35  km  for  Ax  =  300  km,  and  146  km  for  Ax  *  6 00  km, 
when  A0/dp  =  -  6  deg.  (100  mb)  \ 

It  may  be  seen  from  fig.  5f  for  example,  that  the  distance  between 
the  solid  and  the  dashed  curves  does  not  in  itself  indicate  how  6erious 
the  distortion  is.  Therefore,  the  percentage  change  in  the  growth  rate 
of  the  finite-difference  case,  as  compared  to  the  continuous  case,  would 
provide  a  better  measure  of  the  distortion.  Now,  let  5  be  an  arbitrary 
quantity,  and  let  S'  differ  from  ?  by  an  amount  A£.  Then  the  percentage 
change  of  p1  relative  to  £  may  be  written  as 

»  100 • 

This  expression  has  been  used  in  the  construction  of  figs.  10-12.  The 
implication  of  fig.  10  is  that  the  percentage  changes  due  to  the  intro¬ 
duction  of  a  fairly  large  space  mesh  Ax  are  probably  not  serious,  except 
possibly  in  the  immediate  vicinity  of  the  critical  curve.  The  changes 
due  to  a  mesh  size  of  J00  km  amount  to  about  1/3  to  1/2  of  those  shown  in 
the  figure.  Smaller  (larger)  static  stabilities  AG/flp  increase  (reduce) 
those  changes  proportionately.  It  is  well-known  that  effects  due  to  the 
time  mesh  At  are,  in  general,  smaller  than  those  due  to  the  space  mesh  Ax 
when  equations  (1)  and  (2)  are  replaced  by  finite  differences.  This  is 
also  borne  out  by  figs.  11  and  12,  which  indicate  the  relative  insensitivity 
of  the  doubling  times  due  to  the  changes  in  At  alone.  Another  feature  shown 
in  figs.  10-12  is  that  effects  of  Ax  and  At  on  the  growth  rates  are  virtually 
independent  of  the  vertical  shear  U*  over  large  portions  of  the  instability 
region. 
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Figure  11.  Percentage  change  in  doubling  time  due  to  increasing  At  from 
Yz  hr  to  1  hr.  Here  the  latitude  is  45  deg  N,  d0/dp  =  -6  deg  (100  mb)~\ 
Ax  =  300  km,  and  the  units  are  10  ^  per  cent.  The  critical  curve  is  that 
for  Ax  =  300  km. 
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Figure  12.  Same  as  fig.  11 t  but  for  an  increase  in  At  from  ft  hr  to  2  hr. 
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In  contrast  to  the  differential  system,  the  doubling  time  in  the 
finite-difference  case  depends  on  the  undisturbed  zonal  flow  U  due  to  the 
presence  of  in  (16).  This  dependence,  however,  is  slight.  The  per¬ 
centage  changes  in  the  growth  rates  due  to  increasing  U  from  15  to  30  m 
sec~^  are  generally  less  than  ft  per  cent  for  9%/dp  <  0,  and  lees  than 
2  per  cent  in  the  neutral  case,  assuming  Ax  =  300  km,  At  *  1  hr.  Here 
the  effect  of  increasing  Ax  is  to  reduce  the  changes  brought  about  by 
increasing  the  speed  of  the  horizontal  mean  flow. 

4.  Conclusions 


The  introduction  of  implicit  finite  differences  into  the  governing 
equations  of  a  linearized  two-parameter  model  causes  only  slight  distortions 
of  the  growth  rates  of  baroclinically  unstable  waves,  except  in  the  immediate 
vicinity  of  the  critical  curve.  For  the  commonly  used  space  and  time  mesh 
sizes,  Ax  ■  300  km.  At  *  1  hr,  the  percentage  change  in  growth  rate,  com¬ 
pared  to  that  of  the  differential  system,  amounts  to  only  about  2  per  cent 
when  dg/dp  *  -6  deg  (100  mb)~\  and  the  wavelength  L  «  4500  km.  In  general, 
the  distortion  is  greater  for  smaller  static  stabilities,  shorter  waves  and 
larger  Ax.  It  is  very  insensitive  to  changes  in  At,  and  depends  only  slightly 
on  the  undisturbed  zonal  flow.  With  the  exception  of  the  neutrally  stable 
case,  the  growth  rate  of  shorter  waves  is  increased,  and  that  of  the  longer 
waves  slightly  decreased  when  an  implicit  scheme  is  introduced,  as  com¬ 
pared  to  the  continuous  case.  Finally,  the  percentage  changes  in  the 
growth  rates  due  to  changes  in  Ax,  At,  and  U  sure  virtually  independent  of 
the  vertical  shear  U*. 
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appendix 


The  expressions  for  the  real  and  imaginary  parts  of  the  numerical 
phase  speed  C^f  equations  (15)  and  (16),  can  be  obtained  from  (9)  using 
the  relationship 
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(Churchill,  i960),  where  w  =  u  +  iv.  Let 
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and  where  the  upper  sign  corresponds  to  subscript  lf  and  the  lower  sign 
to  subscript  2.  Identity  (18)  now  becomes 
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Here  -  n<  <  *,  and  IntRj/Rg)  is  the  natural  logarithm  of  the  real 

positive  number  Sj/Rg.  Writing  now 
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setting  w  =  XZ,  u  =  XZ^,  v  =  XZ^,  where  X  is  given  by  (8),  and  Z and 
Z^  by  (13)  and  (14)  respectively,  expressions  (15)  and  (16)  follow  directly 
from  (19). 

It  can  also  be  shown  that  C^,  CN  and  CN  tend  to  their  corresponding 

r  1 

counterparts  of  the  continuous  system  when  Ax  and  At  tend  to  zero  simultaneously. 
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